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Summary 


Correction of airfoil data for sidewall boundary-layer 
effects requires a knowledge of the boundary-layer displacement 
thickness and the shape factor with the tunnel empty. To 
facilitate calculation of these quantities under various test 
conditions for the Langley 0.3-m Transonic Cryogenic Tunnel, a 
computer program has been written. This program reads the 
various tunnel parameters and the boundary-layer rake total head 
pressure measurements directly from the Engineering Unit tapes to 
calculate the required sidewall boundary-layer parameters. 

Details of the method along with the results for a sample case 
are presented. 


INTRODUCTION 

The presence of boundary-layers on the sidewalls of a two- 
dimensional wind tunnel are known to affect the test data on the 
airfoils. Recently, methods have been proposed to account for 
this sidewall boundary-layer interference. A summary of these 
correction methods may be found in reference (1) . For applica- 
tion of these methods, it is necessary to have a knowledge of 
the values of the empty test-section sidewall boundary- layer 
displacement thickness and the shape factor. These values are 
generally deduced from the total head pressure measurements 
across the boundary-layer, during the calibration phase of the 
tunnel . 

For the Langley 0.3-m Transonic Cryogenic Tunnel, the side- 
wall boundary-layer measurements were recently made using a wall 
mounted boundary-layer rake in the newly installed adaptive wall 
test-section. To calculate the boundary-layer displacement 
thickness and the shape factor from these measurements, a Fortran 
computer program was written. The present report gives a brief 
description of the program developed and its use for calculating 
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the required sidewall boundary-layer characteristics directly 
from the tunnel data stored on the Engineering Unit tapes. 


Nomenclature 


Test-section width 
Boundary-layer shape factor 
Length scale 
Mach number 

Index in power-law for the velocity profile (Uccy 1 / 11 ) 
Total pressure 
Prandtl number 
Static pressure 

Boundary-layer recovery factor 

Static temperature 

Velocity 

Coordinate normal to the sidewall 


Ratio of specific heats 
Boundary-layer thickness 
Boundary-layer displacement thickness 
Boundary- layer momentum thickness 


Subscripts 


refers to conditions outside the boundary-layer 
refers to boundary-layer 
refers to freestream conditions 


Description of the Method 


The details of calculation of the boundary-layer displace- 
ment thickness and the shape factor are given in reference (2) . 
First, the measured total head pressures in the boundary-layer 
are converted to Mach numbers using the isentropic relation, and 
assuming the static pressure across the boundary-layer to be 
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constant at the freestream value. The isentropic relation for 
the Mach number is given by (for 7 = 1*4) 

j 

M bl - (5{(P bl /p)* 2857 - 1>]°* 5 (1) 

Knowing the Mach number across the boundary- layer, the 
velocity variation can be calculated using the relation 

( U bl/ U e) " (T b l/T e ) 1/2 (M bl /M e ) (2) 

The temperature measurements are generally not made, and 
hence it is necessary to assume a suitable variation acrosss the 
boundary-layer. For the adiabatic wall, the relation between the 
velocity and temperature profiles can be expressed by the 
Crocco's relation (reference 3) given by (for 7= 1.4), 

(Tbi/ T e ) = 1 + 0.2rM e 2 [l - (U bl /U e ) 2 ] (3) 

where r is boundary- layer recovery factor. The value of r 
depends on the nature of the boundary-layer, and the value of 
the Prandtl number (Pr) . For the turbulent boundary-layer, 
this value is approximately Pr 1 / 3 . Using Crocco's relation, 
the velocity and temperature variation across the boundary-layer 
can be calculated. The boundary-layer displacement thickness, 
momentum thickness and the shape factor are then calculated by 
integrating the expressions: 



5 - /[l - ( T e/ T bl>( u bl/ u e)3 <*y 

•l 

(4) 


V /* v 

« - /< T e /T bl )(O bl /U e )[l - <U bl /U e )] dy 

(5) 

and 

H - (5*/0> 

(6) 


For calculating the incompressible values of the displace- 
ment and momentum thickness, the factor (T e /T bl ) is set equal to 
unity, in equations (4) and (5) . 
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Program Description 


The method described above has been incorporated into a 
Fortran computer program called BOUNLAY. The program reads 
the required input data either from the Engineering Unit tapes 
or from user input. The boundary- layer rake tubes are connected 
to scanivalve. The values of the tunnel stagnation pressure, 
stagnation temperature and the static pressure are read at each 
instant the boundary-layer pressure is measured. These instan- 
taneous values are used in the calculation of the boundary-layer 
parameters. 

The boundary-layer rake used in the 0.3-m TCT (figure 1) 
has fifteen tubes and the total height is about 0.6". The tube 
closest to the wall is about 0.024" from the wall. Hence, 
to evaluate the boundary-layer thickness, it is necessary to 
extrapolate the velocity profile up to the wall in a suitable 
manner. Also, under certain flow conditions, it is probable 
that the outermost tube may not be positioned within the 
boundary-layer edge. This requires suitable extrapolation of the 
measured data from the outermost tube to the outer edge of the 
boundary-layer. Hence, a power-law variation was assumed for the 
turbulent boundary- layer velocity profile. The velocity profile 
was represented by the equation 

< u bl/ u e> “ (Y/6) 1/n (6) 

or log(U bl /U e ) = ( 1/n) log (y/L) + (1/n) log(L/6 ) (7) 

where •$ ' is the boundary-layer thickness, 'n' is the index of 
the power-law variation and L is a length scale. The values of 
’6' and 'n* are determined by making a least square fit to the 
calculated values of (U bl /U e ) at different heights across the 
boundary-layer. The least square fit is accomplished by using 
the library subroutine LSQPOL. An example of this least square 
fit result is shown in figure 2. The slope and the intercept 
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of the least square fit determine the values of 'n' and 'S' 
respectively. 

After determining the index of power-law variation for the 
velocity, the integrals for the displacement thickness and the 
momentum thickness are evaluated by using the library subroutine 
CADRE . 

It was recognised that the readings of some of the tubes 
may be erroneous due to extraneous effects during tests. These 
measurement points need to be omitted to avoid erroneous calcula- 
tion of the sidewall boundary-layer parameters. Hence, provision 
has been made to skip up to three tube locations. 

Program Input 

The input data to the program BOUNLAY is through the 
namelist PARAM and the relevant parameters with their default 
values are given in Appendix A. The value of the parameter 
IDATA, either 'O' or ' 1', specifies whether the required tunnel 
parameters are to be read from the Engineering Unit tapes, or 
are given directly as input by the user. The tube heights from 
the wall are specified through the parameter Y. The default 
values shown correspond to the recent measured data, and can be 
changed through the namelist, if new measurements are made, or 
if a new rake is used. Corresponding to the default option of 
IDATA=0 , if the tunnel parameters are to be read from the 
Engineering Unit tapes, the values of the array numbers where the 
various parameters are stored are to be specified through the 
namelist. These parameters and other control parameters required 
are indicated in Appendix A. 

A typical set of input data along with the control cards 
required are shown in Appendix B. Following the namelist PARAM, 
the test, run and point numbers are specified in 313 format, when 
the tunnel parameters are to be read from the Engineering Unit 
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tapes. Specifying 'O', for the point number allows the program 
to calculate the sidewall boundary- layer parameters for all the 
points in the specified run number. For termination of the cal- 
culation, the run number is specified as '999', in the last 
record . 

If the tunnel parameters are given directly as input data, 
following the namelist PARAM, the program reads a title card in 
the format 10A8. After this, the values of the stagnation 
pressure (psia) , stagnation temperature (degree K) , height of the 
tube from the wall (inches) , boundary-layer pressures (psia) and 
tunnel static pressure (psia) are specified in format 5F10.6. 
These are repeated corresponding to the value of NBL specified 
in the namelist PARAM. 

Results for a Sample Case 

The use of the computer program developed is demonstrated 
by applying it to one of the runs conducted recently during empty 
test section calibration (Test 200) . The sample output shown in 
Appendix C, corresponds to the results obtained with the data in 
Appendix B. The corresponding nominal freestream Mach number 
and Reynolds number are 0.7 and 30 million/ft, respectively (Run 
number 19 and Point number 6) . 

First, the values of the various parameters in namelist 
PARAM are printed. Following this, the title, and the test, 
run and point numbers are printed. Subsequently, the tunnel 
conditions averaged over the time interval of the boundary- 
layer measurements, the calculated values of the boundary- 
layer parameters are printed. At the end, the variation of 
the pressure, temperature, velocity and Mach number across the 
boundary-layer are tabulated. 

If the value of the parameter IPLOT is not equal to zero, 
the calculated velocity profiles are written on a temporary file 
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on Tape 10. These are then plotted using the program UPLOT, and 
the corresponding control cards required for plotting are shown 
in Appendix A. The program UPLOT uses the library subroutine 
INFOPLT and plots both the measured velocities and the calculated 
power-law profile for comparison (see Figure 3) . 

For the sample case considered, it is seen that the value of 
the displacement thickness and the shape factor at the model 
location are .0713" and 1.4595 respectively. For the testsection 
width of 13", the value of the parameter (2 5*/k>) which is one of 
the parameters characterizing the extent of sidewall boundary- 
layer correction is about 0.01. The corresponding corrrection 
for the test Mach number on long chord models is about -0.008. 
However, for short chord models, it is likely the sidewall 
boundary-layer correction will be much smaller. The detailed 
calculation of the sidewall boundary-layer parameters, and the 
corresponding corrections over the operating range of the 0.3-m 
TCT will be dealt in a future report. 
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I DATA 
I COMP 
I PLOT 
L 

NBL 

Y 

ISKIP 

IBL 

I PRINT 
IOPT 


NTEST 

NRUN 

NPOINT 

NPORT 

NBLSCAN 

NPTINF 

NPRT 

NPTC 


Appendix A 

Parameters in Namelist PARAM 

= 0, Tunnel parameters are read from EU tapes, (0)* 

= i. Tunnel parameters are given as input data 
= 0, Incompressible values of 5* and Q calculated, (1) 
= l, Compressible values of 6 and 6 calculated 
= 0, Plot files are not generated, (1) 

= 1, Plot files are generated 
Reference length scale, (1.0) 

Number of tubes in the boundary- layer rake, 

Maximum value =30, (15) 

Tube heights from the wall (.024, .06, .10, .14, .18, 

.22, .26, .30, .34, .38, .42, .46, .50, .55, .59) 

Number of tubes to be omitted calculation, Max. 3 (1) 
Port number of the tubes to be omitted (1, 0, 0) 

= 0, Only final results are printed (0) 

= 1, Intermediate results also printed 
= 1, 6* and 6 calculated using the power-law for 
the entire profile. 

= 2 , 5* and 6 calculated using the power-law only 

upto first tube, and using trapezoidal rule 
for the remaining profile 
Array number for the test + (3) 

Array number for the run + (5) 

Array number for the point 4- (6) 

Array number for the port + (10) 

Array number for the boundary-layer 
rake pressure + (117) 

Array number for the tunnel total pressure* (107) 
Array number for the tinnel total temperature + (217) 
Array number for the tunnel static pressure* (108) 


* Numbers within paranthesis refer to default values 
+ to be specified only when IDATA = 0. 
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Appendix B 

Control Cards and Sample Input Data 


PROBOUN,T1000, CM170000. Bldgnum name 

USER , usernum . passwor . 

CHARGE , chgnum , LRC . 

ATTACH , TAPE7=D* ***** */UN=usernum . 

GET , BOUNLAY . 

FTN , I = BOUNLAY , L=0 , PL=10000 . 

ATTACH ( FTNMLIB/UN=LIBRARY , NA) 

LDSET ( LI B=FTNMLI B , PRESETA=INDEF) 

HAP, OFF. 

LGO. 

REWIND, TAPE10. 

RETURN (*,TAPE10) 

GET , UPLOT 

FTN , I=UPLOT , L=0 . 

ATTACH , LRCGOSF/UN-LIBRARY , NA . 

LDSET ( LI B=LRCGOS F , PRESETA=INDEF) 

LGO. 

DELIVER. bldgnum name 
PLOT . VARIAN , FF 
/eor 

$PARAM IDATA=0 , IOPT=2 , $END 
200019000 
200999000 

$PARAM IDATA=1 , IOPT=2 , $END 


TEST NO: 

200. RUN: 

19. 

POINT: 

6. SAMPLE INPUT DATA 

35.5818 

129.6503 

.0240 

29.1408 

25.6570 

35.5475 

129.8012 

.0600 

30.9015 

25.6787 

35.5341 

129.8498 

.1000 

31.5694 

25.5685 

35.6083 

129.8417 

.1400 

31.9734 

25.4811 

35.6049 

129.8197 

.1800 

32.5464 

25.5275 


continued 
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35.6106 

129.7544 

.2200 

32.8914 

25.5226 

35.6429 

129.6369 

.2600 

33.5220 

25.5598 

35.6908 

129.5605 

.3000 

33.8837 

25.6194 

35.5991 

129.6841 

.3400 

34.2634 

25.6547 

35.5794 

129.7295 

.3800 

34.8155 

25.5787 

35.5810 

129.7403 

.4200 

34.8841 

25.5163 

35.6308 

129.7172 

.4600 

35.1331 

25.5355 

35.5785 

129.7094 

.5000 

35.2827 

25.5193 

35.6037 

129.6764 

.5500 

35.4686 

25.5294 

35.6374 

129.5828 

.5900 

35.4827 

25.4787 


/eor 

/eof 
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Appendix C 

Sample Output from the Program 


$PARAM 




NTEST 

= 

3 , 


NRUN 

= 

5, 


NPOINT 

= 

6, 


NBL 

= 

15, 


NPORT 

= 

10, 


NBLSCAN 

= 

117, 


NPTINF 

= 

107, 


NPTC 

= 

108, 


NPRT 

= 

217, 


ICOMP 

= 

1, 


I PLOT 

= 

1, 


Y 

= 

. 24E-01, 

. 6E-01 , 



.22E+00, 

« 26E+00 



. 42E+00 , 

. 46E+00 



H 

* 

H 

H 

H 

H 

H 

L 

= 

. 1E+01 , 


ISKIP 

= 

1, 


IBL 

= 

1, 0, 0, 


I DATA 

= 

0, 


I PRINT 

= 

1, 


IOPT 

S 

2, 


$END 





.1E+00, . 14E+00 , . 18E+00 , 

. 3E+00, .34E+00, .38E+00, 

.5E+00, . 55E+00 , .59E+00, 
I, I, I, I, I, I, I, I, I, 


continued 
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Sample Output from the Program 


0.3-M TCT STREAMLINED WALL TEST SECTION 
RIGHT SIDEWALL BOUNDARY-LAYER MEASUREMENTS 


TEST: 200 RUN: 19 POINT: 6 


AVERAGE MACH NUMBER 
AVERAGE REYNOLDS NUMBER/ FT 
AVERAGE STAGNATION PRESSURE PSIA 
AVERAGE STAGNATION TEMPERATURE K 
PRANDTL NUMBER 
RATIO OF SPECIFIC HEATS 
RECOVERY FACTOR 
REFERENCE LENGTH SCALE L 


.7049 

.295E+08 

35.6035 

129.7217 

.7600 

1.4000 

.9126 

1.0000 


VALUE OF N IN POWER LAW PROFILE : 7.9196 
BOUNDARY-LAYER THICKNESS DELTA/L : .5642 

DISPLACEMENT THICKNESS DELTASTAR/L : .0713 
MOMENTUM THICKNESS THETA/L : .0488 
SHAPE FACTOR : 1.4595 


A: DISPLACEMENT AND MOMENTUM THICKNESS REFER 
TO COMPRESSIBLE VALUES. 

B: DISPLACEMENT AND MOMENTUM THICKNESS CALCULATED 
USING POWER-LAW FROM WALL TO FIRST TUBE. 


continued 
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NO. 

Y/L 

M/MINF 

1 

.0600 

.7469 

2 

.1000 

.7935 

3 

.1400 

.8172 

4 

.1800 

.8489 

5 

.2200 

.8677 

6 

.2600 

.8991 

7 

.3000 

.9149 

8 

.3400 

.9372 

9 

.3800 

.9650 

10 

.4200 

.9684 

11 

.4600 

.9777 

12 

.5000 

.9868 

13 

.5500 

.9940 

14 

.5900 

.9932 


P/PT 

T/TINF 

U/UINF 

8693 

1.0374 

.7608 

8884 

1.0315 

.8059 

8979 

1.0287 

.8288 

9141 

1.0239 

.8590 

9236 

1.0211 

.8767 

9405 

1.0162 

.9063 

9494 

1.0137 

.9211 

9625 

1.0101 

.9419 

9785 

1.0057 

.9678 

9804 

1.0052 

.9709 

9860 

1.0037 

.9795 

9917 

1.0022 

.9878 

9962 

1.0010 

.9945 

9957 

1.0011 

.9938 
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Appendix D: Program listing 


PROGRAM CRYOTBL(INPUT, OUTPUT, TAPE5=INPUT,TAPE7,TAPE10) 

c • *•••*••••*••*»•**•**•***** 

C 0.3-M TCT BOUNDARY • LAYER CALCULATION 

C 

C THIS PROGRAM CALCULATES THE SIDEWALL BOU N D AR Y - L A Y E R 

C DISPLACEMENT THICKNESS, MOMENTUM THICKNESS AND THE 

C SHAPE PARAMETER FROM THE TOTAL HEAD PRESSURE MEASURE- 

C MENTS IN THE BOUND ARY • L AYER . 


C 

REAL HBL,MBLBMIN,MINF,L,MAV 
DIMENSION XX(30),YY(30),EPS(3),YM(30) 

DIMENSION U(30) , RE S I D ( 30 , 1 ) , SUM ( 1 ) , A( 30 , 1 ) , B ( 3 0 , 1 ),UK(30,2) 
DIMENSION Y(30) ,PBL(30) , PS(30) , ARRAY ( 1000) 

DIMENSION P1(30),P2(30),P3(30),P4(30),P5(30) 

DIMENSION PTINF(30),TTINF(30),MBL(30) 

DIMENSION PBLNON<30),MBLBMIN(30),UBLBUIN(30>,TBLBTIN<30) 
DIMENSION U(30),T (30), F(30),G(30), I B L ( 3 ) 

DIMENSION P(30),YN(30),PN(30),M(30) 

DIMENSION IORF1(30),PTC(30),PRT(30) 

DIMENSION HEAD I NG ( 1 0 > 

NAMELIST /PARAM/ N T E S T , N R UN , N PO 1 N T , N B L , N POR T , N B L S C A N , N P T I N F , 
SNPTC.NPRT, ICOMP, IPLOT,Y,L, ISKIP, IBL, IDATA, IPRINT , IOPT 

600 CONTINUE 

C SET DEFAULT VALUES OF THE NAMELIST PARAMETERS 


IDATA 

= 

0 

NTEST 

= 

003 

N RUN 

= 

005 

NPOI NT 

= 

006 

NBL 

= 

015 

NPORT 

= 

010 

NBLSCAN 

s 

117 

NPRT 

= 

217 

NPT I NF 

= 

107 

NPTC 

= 

108 

I COMP 

= 

001 

I PLOT 

= 

001 

Y< 1 ) 

= 

0.024 

Y( 2) 

= 

0.060 

Y ( 3) 

= 

0.100 

Y ( 4) 

s 

0.140 

Y ( 5) 

= 

0.180 

Y ( 6) 

s 

0.220 

Y ( 7) 

s 

0.260 

Y ( 8) 

= 

0.300 

Y ( 9) 

= 

0.340 
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Y( 10) 

Y<11> 

Y( 12) 

Y ( 1 3 ) 

Y( 14) 

Y < 1 5 ) 

L 

I SKI P 
IPRINT 
I B L ( 1 ) 

I B L ( 2 ) 

I B L ( 3 ) 

I OPT 

C SET ERROR FLAG TO ZERO 

IERR = 0 


C READ INPUT DATA FROM NAMELIST 

140 R E AD ( 5 , PARAM) 

IF(EOF<5)) 120, 130 

120 STOP 120 
130 CONTINUE 

I F ( I D A T A .EG. 0) PRINT PARAM 

C STORE Y VALUES 

DO 1600 II = 1 , N B L 
YM ( I I ) = Y ( I I ) 

1600 CONTINUE 

ENBL = NBL 
ISKP = ISKIP 


I F ( IDATA .EQ. 0) GO TO 700 

C READ TITLE CARD 

READ 650, <I<EADING< I ), 1-1 ,8) 

C READ INPUT DATA 

DO 710 1=1 , NBL 

READ 720, PT I N F< I ) , TT INF ( I ), Y< I ) , PBL< I ) ,PTC< I ) 

PRINT 720 , P T I N F ( I ),TTINF(I),Y(I ) , P B L ( I ),PTC(I ) 

720 FORMAT (5F10.6) 

710 CONTINUE 

PRINT 660, (HEADJNG(I), 1-1,8) 

GO TO 15 

C TUNNEL AND B OU N D A R Y - L A Y E R PARAMETERS ARE READ FROM EU TAPES 

700 CONTINUE 


* 0.380 
= 0.420 
= 0.460 
= 0.500 
= 0.550 
= 0.590 
= 1.000 
= 1 
= 0 
= 1 
= 0 
= 0 
= 1 
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c FIRST TWO FRAMES CONTAIN NON-DAT, INFOR SUCH AS UNITS AND LABELS 

DO 1000 l 1«1 ,2 

READ(7)KEY,N, (ARRAY( I ), 1*1 ,N) 

1000 CONTINUE 

C SET REQUIRED INPUT DATA TO VARIABLES USED IN THE PROGRAM 

N S V 0 = NBLSCAN 


2 READ 4, IT, IR, IP 

I F ( I R . EQ . 999 ) GO TO 110 

C READ THE TUNNEL VARIABLES FROM THE EU TAPE 

5 READ(7)KEY,N, (ARRAY( I ), 1*1, N) 

I F(EOF(7) )2,6 

6 CONTINUE 

C SELECT THE REQUIRED PARAMETERS 

IF (IP. NE . 0) GO TO 1500 
I F ( ARRAY ( NRUN ) .GT. IR) GO TO 2 

1500 CONTINUE 

1 TEST*ARRAY ( NTEST ) 

!RUN=ARRAY(NRUN) 

I POI NT*ARRAY ( NPOI NT ) 

I F ( I RUN .EQ. 0) GO TO 7 
I F( I RUN . LT . I R ) GO TO 5 
I F ( I P .EQ . 0 ) GO TO 7 
I F ( I POINT .NE. IP) GO TO 5 

7 I PORT = A R R A Y (NPORT ) 

I F ( I PORT .EQ. 0) GO TO 5 
I F ( I PORT . GE . 16) GO TO 5 

IORF1(IPORT)=ARRAY(NPORT) 

P B L ( IPORT )=ARRAY(NSVO) 
PTINF(IPORT)=ARRAY(NPTINF) 

P T C ( IPORT ) = ARRAY(NPTC) 
TTINF(IPORT)=ARRAY(NPRT) 

IF ( I PORT .EQ. 1 ) GO TO 10 
I F ( IPORT .EQ.NBL) GO TO 15 
GO TO 5 

10 PRINT 40 

PRINT 100, ITEST, IRUN, IPOINT 

13 GO TO 5 
15 CONTINUE 

I F( I SKI P .NE . 0) 
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SC ALL SKIP(NBL, ! SK I P , IBL,Y,PBL,PTINF,TTINF,PTC, IPRINT , IERR) 
S I GMA = 0 . 76 
GAMMA* 1 .4 
N = NBL 

C CALCULATE B OU NO A R Y • L A Y E R PROFILES 

CALL BLPROFL(NBL, Y.PBL, PTC, PT INF, TT INF, L, GAMMA, SIGMA, 
SYM,PBLN0N,MBLBHIN,UBLBUIM,TBLBT1N,MAV,PTINFAV,TT1NFAV, 
SREYAV.REC, IPRINT, IERR) 

I F ( I E R R .ME. 0) GO TO 600 

C FIT A POWER-LAW CURVE TO THE VELOCITY PROFILE 

CALL LSFIT(NBL,YN,UBLBUIN,L, IPRINT, EN, DELTA, IERR) 

I F ( IERR .NE . 0) GO TO 600 

C CALCULATE DISPLACEMENT THICKNESS AND MOMENTUM THICKNESS 

CALL THIKNES(NBL,YN,UBLBU!N,TBLBTIN,L, ICOMP, DELTA, EN, IOPT 
S,MAV,REC, GAMMA, DELS, THETA, IPRINT, IERR) 


PRINT 

FINAL 

RESULTS 

PR I NT 

500, 

MAV 

PRINT 

501 , 

RE YAV/3 . 2808 

PRINT 

502, 

PT I N FAV 

PRINT 

503, 

TTINFAV 

PRINT 

504 , 

S I GMA 

PRINT 

505 , 

GAMMA 

PRINT 

506, 

REC 

PRINT 

507, 

L 

PRINT 

508, 

EN 

PRINT 

509, 

OELTA/L 

PRINT 

510, 

DELS 

PRINT 

511, 

THETA 

PRINT 

512, 

DELS/THETA 


I F ( I COMP . EQ . 

0) 

PRINT 

521 

I F ( I COMP . NE . 

0) 

PRINT 

522 

IFCIOPT .EQ. 

1 ) 

PR I NT 

523 

I F ( IOPT .EQ. 

2) 

PRINT 

524 


PRINT 530 

DO 550 I = 1, NBL 

PRINT 531, I , YN< I ) ,MBLBMIN< I ) .PBLNONC I ),TBLBT INC I ) ,UBLBUIN< I ) 
550 CONTINUE 

C WRITE PLOT FILE 

I F ( I P L OT .EO. 0) GO TO 570 
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WRITEdO, 533) NBL, EN 
I F ( I D A T A .EQ. 0) WRITE(10 
I F ( I DATA .NE. 0) WRITEdO 


534) ITEST, I RUN , I PO I NT 

535) HEADINGd >,HEADING(2) 


DO 560 I » 1 , NBL 

WRITEdO, 536) Y(I) /DELTA, UBLBU I N ( I ) 
560 CONTINUE 


570 PRINT 532 

I F ( I DATA .NE. 0) GO TO 600 
NBL - ENBL 
ISKIP = ISKP 
C RESET Y VALUES 

DO 1 700 It > 1 , N B L 
Y ( I I ) = Y M C I I ) 

1700 CONTINUE 
GO TO 5 


500 FORMAT !/10X, ‘AVERAGE MACH NUMBER 

501 FORMAT ( 1 OX , * AVERAGE REYNOLDS NUMBER/FT 

502 FORMAT ( 10X, ‘AVERAGE STAGNATION PRESSURE PSIA 

503 FORMAT ( 10X, ‘AVERAGE STAGNATION TEMPERATURE K 

504 FORMAT! 10X,*PRANDTL NUMBER 

505 FORMAT! 10X, ‘RATIO OF SPECIFIC HEATS 

506 FORMAT! 1 OX , * R E COVE R Y FACTOR 

507 FORMAT! 1 0 X , * R E F E R E N C E LENGTH SCALE L 

508 FORMAT!/10X, ‘VALUE OF N IN POWER LAW PROFILE 

509 FORMAT! 1 OX , ‘BOUNDARY - LAYER THICKNESS DELTA/L 

510 FORMAT!/10X, ‘DISPLACEMENT THICKNESS DELTASTAR/L 

511 FORMAT! 1 0 X , * M OM E N T U M THICKNESS THETA/L 

512 FORMAT! 10X, ‘SHAPE FACTOR 


. * 

F10.4) 

. * 

E10.3) 

. * 

F10.4) 

. * 

F10.4) 

• ♦ 

F10.4) 

• * 

F10 .4) 

- * 

F10.4) 

. * 

F10.4) 

. * 

F10.4) 

. * 

F10.4) 

. * 

F10.4) 

: * 

F10.4) 

. * 

F10.4) 


521 FORMAT !//1 OX , *A : DISPLACEMENT AND MOMENTUM THICKNESS REFER*, 

S/10X, • TO INCOMPRESSIBLE VALUES.*) 

522 FORMAT(//10X,*A: DISPLACEMENT AND MOMENTUM THICKNESS REFER*, 

S/10X, * TO COMPRESSIBLE VALUES.*) 

523 F ORMAT ! 1 OX , *B : DISPLACEMENT AND MOMENTUM THICKNESS CALCULATED*, 

$ /10X,* USING POWER-LAW FOR WHOLE PROFILE.*) 

524 FORMAT ! 10X , ‘B : DISPLACEMENT AND MOMENTUM THICKNESS CALCULATED*, 

$ /10X,* USING POWER-LAW FROM WALL TO FIRST TUBE.*) 


530 FORMAT!///10X,*NO.*,6X,*Y/L*,4X,*M/MINF*,6X,*P/PT*,4X, 
S*T/TINF*,3X,*U/UINF*,/) 

531 FORMAT(10X, I2.5F10.4) 

532 FORMAT! //30X, * *) 

533 FORMAT!I5,F10.6) 

534 FORMAT !*17HU/UE *,314,* Y/DELTA *) 

535 FORMAT < * 1 7 H U/UE *,2A10,* Y/DELTA *) 

536 FORMAT!2F10.6) 
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650 FORMAT (8A10) 

660 FORMAT(1H,/10X,8A10,/> 

4 FORMAT(3I3) 

40 FORMAT(1H1/10X, *0.3-M TCT STREAMLINED WALL TEST SECTION*, 
S/10X, ‘RIGHT SIDEWALL BOUN D AR Y • L A Y E R MEASUREMENTS*) 

100 FORMAT(/10X,*TEST : * , I 3 , 5 X , *RUN : * , I 3 , 5X , * PO I N T : * , I 3 ) 

101 F0RMAT(3I10) 

110 CONTINUE 

GO TO 140 

STOP 

END 
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SUBROUTINE SKIP(NBL, ISKIP, IBL, Y,PBL,PTINF,TTINF, PTC, IPRINT.IERR) 

DIMENSION IBL(3),Y(30),PBL(30),PTINF(30),TTINF(30),PTC(30) 
DIMENSION P1(30),P2(30),P3(30),P4(30),P5<30) 

Q ********************** 

C 

C NBL : NUMBER OF TUBES IN BOUNDARY - LAYE R RAKE 

C ISKIP : NUMBER OF POINTS TO BE SKIPPED (MAXIMUM 3) 

C I B L ( 3 ) : PORT NUMBERS TO BE SKIPPED 

C Y : DISTANCE OF TUBES FROM WALL 

C PBL : TOTAL PRESSURE IN THE BOUN D A R Y • l A Y E R 

C PTINF : FREE-STREAM STAGNATION PRESSURE 

C TTINF : FREE-STREAM STAGNATION TEMPERATURE 

C PTC : FREE-STREAM STATIC PRESSURE 

C IPRINT : NOT USED 

C I ERR : NOT USED 

C 

C NOTE : ON RETURN, NEW VALUES OF N B L , Y , P B L , P T I N F , T T I N F , P T C 

C WILL BE AVAILABLE. 

C 

c ********************** 

JS * 0 

300 DO 100 1*1 , N B L 

DO 200 J = 1 , ISKIP 
I F( I .EO. I B L ( J ) ) GO TO 150 
200 CONTINUE 
GO TO 250 

150 JS * JS ♦ 1 
GO TO 100 


250 II* I • JS 

PI ( I 1 >* P B L ( I ) 
P2( ID* PTINFC I ) 
P3< ID* T T I N F ( I ) 
P4( I D* P T C ( I ) 
P5( I D* Y< I) 

100 CONTINUE 

NBL* NBL - ISKIP 

DO 400 1*1, NBL 
P B L ( I > = P 1 ( I > 

PT I NF ( I) = P2( I ) 
TT I NF ( I ) * P3( I ) 
P T C ( I ) = P4( I) 

Y ( I ) = P5 ( I ) 

400 CONTINUE 


1F(I PRINT .EQ. 0) RETURN 
PRINT 500 
DO 550 1 = 1 , NBL 

PRINT 600, l,Y(I),PTINF(I),TTINF(I),PTC(I),PBL(I) 

550 CONTINUE 

500 FORMAT (/ 1 0X , *RE ARRANGED VALUES OF THE PARAMETERS Y,PTINF,*, 
%/ t *TTINF,PTC,PBL BY SUBROUTINE SKIP ARE AS FOLLOWS*,/) 

600 FORMAT(I10,5F12.6) 

RETURN 

END 
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SUBROUTINE LSFIT(N,Y,U,L,IPRINT # EN, DELTA, IERR) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


******* 


THIS SUBROUTINE FIT A POWER-LAW CURVE TO THE VELOCITY PROFILE TO 
DETERMINE THE CONSTANTS EN AND DELTA 


N 

Y(N ) 

U ( N ) 

L 

I PR I NT 
IERR 


NUMBER OF POINTS IN VELOCITY PROFILE 
VALUES OF Y 

VALUES OF NON-DIMENSIONALIZED VELOCITY 
LENGTH SCALE 

NOT, EQUAL TO ZERO, PRINTS OUTPUT PARAMETERS 

EQUAL TO ZERO, NORMAL RETURN 

NOT EQUAL TO ZERO, ERROR IN CURVE FITTING 


EN : INDEX IN POWER-LAW VELOCITY PROFILE 

C2 : CONSTANT TERM DELTA « L/EXP(C2/C1) 


THE CURVE FIT IS OF THE TYPE 

LOGN ( U/U I N F ) = Cl * LOGN(Y/L) + C2 


********************** 


REAL L 

DIMENSION Y(30),U(30) 

DIMENSION XX(30),YY{30) 

DIMENSION W(30),RES1D(30,1),$UM(1),A(30,1),B(30,1),WK(30,2) 


DO 1 00 I = 1 , N 
XX ( I ) = ALOG( Y( I ) ) 
Y Y ( I ) * A L OG ( U ( I ) ) 
W ( I ) =1 

100 CONTINUE 


C SET PARAMETERS FOR LEAST SQUARE FIT 

NMAX = 30 
N - N 

LL = 1 
MMAX = 4 
M -2 

C LEAST SQUARE FIT FOR N READINGS 

CALL LSQPOL(NMAX,N,XX,LL,YY,W,MHAX,M, R E S I D , SUM , A , B , WK , I E R R ) 

I F ( I ERR .NE. 0) PRINT 500 

500 FORMATC/IOX, ‘ERROR IN SUBROUTINE LSQPOL: CHECK VELOCITY DATA*) 

I F( IERR .NE. 0) RETURN 
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C DETERMINE EN AND DELTA 

EN = 1 ./B(2,1 ) 

DELTA * L/EXP(B(1,1)/B(2,1)) 

IF < IPRINT .EO. 0) GO TO 300 

C PRINT INPUT TO AND OUTPUT FROM SUBROUTINE LSQPOL 

PRINT 1 000, NMAX , N , L L , MMAX , M 
DO 200 1= 1, N 

PRINT 2000, I,W(I),Y(I),U(I),XX(I),YY(1),RE$ID(I,1) 

200 CONTINUE 

PRINT 2020 , I ERR , EN, DELTA, SUM(1) 

300 CONTINUE 

1000 F ORMAT ( 1 OX , 5 I 5 ) 

2000 FORMAT < 5 X , l 5 , 6 F 9 . 4 > 

2020 FORMAT(10X,*IERR =*,I2,5X,*EN ■ * , F8 . 4 , 5X , *D E LT A«* , F 8 . 4 , 
S5X,*SUM(1 ) = *, F8.4) 

2030 FORMAT(/10X,F10.6,/10X,5E10.4> 

RETURN 

END 
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SUBROUTINE BLPRO F L ( N , Y , PBL , PS f PT I N F , TT I N F , L , GAMMA , S 1 GMA , 
$YN,PBLNON,MBLBMIN,UBLBUIN,TBLBTIN,MAV,PTINFAV,TTINFAV, 
$REYAV,REC, IPRINT, IERR) 

REAL MBL,MBLBMIN,M1NF, L , MA V 
DIMENSION Y<30) ,PBL<30) ,PS<30) 

DIMENSION PTINF<30),TTINF(30),MINF(30),MBU30) 

DIMENSION YN<3Q),PBLNON(30>,MBLBMIN<30) f UBLBUIN(30),TBLBTIN<30> 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


********************* 

N : NUMBER OF TUBES IN THE BOUNDARY -LAYER RAKE 

Y s PROBE HEIGHTS FROM THE WALL 

PBL : TOTAL PRESSURE IN THE BOUND AR Y * LA YER 

PS : STATIC PRESSURE ( F R E E - S T R E AM ) 

PTINF : TUNNEL STAGNATION PRESSURE 

TTINF : TUNNEL STAGNATION TEMPERATURE 

L : LENGTH SCALE FOR NON - D 1 ME N S I ON A L I S I N G Y 

GAMMA : RATIO OF SPECIFIC HEATS 

SIGMA : PRANDTL NUMBER 

YN : VALUES OF Y NORMALISED BY L 

PBLNON : VALUES OF PBL NORMALISED BY PTINF 

MBLBMIN: RATIO OF MACH NUMBER IN BOUND RY - LAYER TO FREE-STREAM 

UBLBU1N: RATIO OF VELOCITY IN BOUNDARY* LAYER TO FREE * STREAM 

TBLBTIN: RATIO OF TEMPERATURE IN BOUND AR Y - L A YE R TO FREE-STREAM 

MAV : AVERAGE MACH NUMBER 

PTINFAV: AVERAGE STAGNATION PRESSURE 

TTINFAV: AVERAGE STAGNATION TEMPERATURE 

REYAV : AVERAGE REYNOLDS NUMBER/METER 

REC : RECOVERY FACTOR 

IPRINT : 

IERR : 

********************** 
RECOVERY FACTOR 

R = SIGMA**(1 . / 3 . ) 

REC = S I GMA** ( 1 . /3 . ) 

CONSTANTS 

G = GAMMA 
61 s (G- 1 . ) / G 

G 2 ® 1./G 
G3= R*(G*1.)/2. 
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C CALCULATE MACH NUMBERS IN BOUND A R T - L A T E R AND FREESTREAM, AND RATIO 

DO 100 I = 1, N 

XI = <PBl< I >/PS< I ) >**G1 -1. 

X2 *(PTINF( I )/PS(! ) ) **G 1 -1. 

I F ( X 1 .LT. 0. .OR. X2 .LT. 0.) GO TO 1000 
XI * <2./<G-1.))*X1 
X 2 * ( 2 . / ( G • 1 . )>*X2 
MB L ( I ) = SORT ( X 1 ) 

M I N F ( I ) = SORT (X2) 

MBLBMIN(I) = MBL( I ) / M I NF< I ) 

PBLNON(I) = PBL(I)/PTINF(I) 

100 CONTINUE 

C CALCULATE RATIO OF TEMPERATURE IN BOUN D AR Y - L A Y E R TO FREESTREAM 

DO 200 I = 1, N 

TBLBTIN(I)= <1. + G3* ( M I N F ( I ) * * 2 . ) > / (1. ♦ G3* ( MB L ( I ) * * 2 . ) ) 

200 CONTINUE 

C CALCULATE RATIO OF VELOCITY IN BOU ND AR Y - L A Y E R TO FREESTREAM 

DO 300 I = 1, N 

UBLBUIN(I) = MBLBMIN(I) • SQ RT ( T B L B T I N ( I ) ) 

300 CONTINUE 

C NON • D I MENS I OAL I ZE RAKE PROBE HEIGHTS TO LENGTH SCALE L 

DO 400 I = 1 , N 
Y N ( I ) = V C I ) / L 
400 CONTINUE 

C CALCULATE AVERAGE QUANTITIES 

SUMM = 0. 

SUMP = 0. 

SUM T = 0. 

DO 500 I = 1, N 

SUMM = SUMM + MINF(I) 

SUMP = SUMP + PTINFd ) 

SUMT = SUMT + TTINF(I) 

500 CONTINUE 

MAV = SUMM/N 

PTINFAV = SUMP/N 
TTINFAV = SUMT/N 

I F( I PR I NT .EQ. 0) GO TO 1010 
PRINT 1020 
DO 1030 I »1, N 

PRINT 1 040, I, MI NF ( I ) , MBL ( I ) ,UBLBUIN( I ) , TBLBT I N( I ) 


25 



1030 CONTINUE 

1 020 FORMAT (/IX,* I MINF MB L U/UINF T/TINF*) 

1040 FORMAT(1X, I5.4F10.6) 

C CALCULATE AVERAGE REYNOLDS NUMBER/METER 

1010 CONTINUE 

C = 1. + <<G - 1 . >/2 . ) • ( M A V* * 2 . ) 

A 1 = PTINFAV/(TTINFAV**1 .4) 

A2 = MAV/< C**2 . 1 ) 

A3 = <64314. /14. 696>*<10. **6. ) 

2 = COMP(PTINFAV,TTINFAV) 

R E Y A V= A3*A2*A1*Z 

I F ( I PR I N T .NE. 0) PRINT 1 050, M A V , P T I N F A V , T T I N F A V . Z , RE Y A V 
1050 FORMAT (/10X,*MAV, PTINFAV, TTINFAV, Z, REYAV: *, /10X, 

S4F10.4.E14.4) 


RETURN 

1000 PRINT 1100 

1100 FORMAT(/10X,*ERROR IN CALCULATING BOUND ARY - LAYER VELOCITY*, 
S* PROFILE: CHECK INPUT PRESSURE DATA*) 

RETURN 

END 
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SUBROUTINE THJKNE$<N,Y,U,T, L, IC OMP, DELTA, EN, IOPT 
$,MINF , RE C, GAMMA, DELS, THETA, I PR I NT, I ERR) 

REAL MINF, L 

DIMENSION Y(30),U(30),T(30), F(30),G(30) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


N 

Y 

u 

T 

L 

I COMP 

DELTA 

EN 

MINF 
REC 
GAMMA 
I PR INT 
IOPT 

IOPT 


I ERR : 


NUMBER OF POINTS IN THE BOUND A R Y - L A Y E R 
RAKE PROBE HEIGHTS 

NON-DIMENSIONALISED VELOCITY PROFILE AT N POINTS 
NON-DIMENSIONALI SED TEMPERATURE PROFILE AT N POINTS 
LENGTH SCALE 

= 0, CALCULATES INCOMPRESSIBLE VALUES 
/ 0, CALCULATES COMPRESSIBLE VALUES 
BOUNDAR - LAYER THICKNESS 

CONSTANT IN POWER-LAW VELOCITY PROFILE 

MACH NUMBER 

RECOVERY FACTOR 

RATIO OF SPECIFIC HEATS 

NOT = ZERO, PRINTS INTERMEDIATE RESULTS 
= 1, DELTASTAR CALCULATED USING EN AND DELTA 
FOR THE ENTIRE PROFILE 

= 2, DELTASTAR CALCULATED USING EN AND DELTA 

ONLY UPTO FIRST TUBE, AND USING TRAPEZOIDAL 
RULE FOR OTHER POINTS. 

NOT USED 


C 

C 


DELS : DISPLACEMENT THICKNESS / L 

THETA : MOMENTUM THICKNESS / L 


C 

c * * 


NBL = N 

I F ( IOPT .EG. 1 ) GO TO 680 


C CALCULATE THE V E L OC I T Y / M A S S DEFECT 

DO 100 1= 1, N 

I F ( I COMP .NE. 0) F ( I ) = 1. - U(I)/T(I) 

IFCICOMP .NE. 0) G ( I ) = ( U ( I ) / T ( I ) ) * ( 1 . - U(I)> 

I F ( I COMP .EO. 0) F ( I ) = 1. - U( I ) 

I F ( I COMP -EG. 0) G < I > = U(I)*(1. - U(I>) 

100 CONTINUE 

C INTEGRATE FROM TUBE(1) TO TUB E ( N ) 

IFCIPRINT .NE. 0 .AND. ICOMP .NE. 0) PRINT 500 
IFdPRINT .NE. 0 .AND. ICOMP .EQ. 0) PRINT 510 
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500 F ORMAT (/IX,* I Y U T (1 • U/T) SUM*, 

$* (U/TX1-U) SUM*) 

510 FORMAT (/IX,* I Y U T (1-U) SUM*, 

$* U(1-U> SUM*) 


SUMI « 0. 

SUMG « 0. 

DO 200 1= 2, N 

1 F ( I .EQ. 2 .AND. IPRINT. NE . 0) 

SPRINT 600, (1-1), Y(1),U(1),T(1),F(1), SUMI, G(I), SUMG 

DY « Y(I) - Y< I - 1 > 

SUM I = SUMI ♦ ( F( I ) + F( I - 1 ) )*(DY/2. ) 

SUMG= SUMG ♦ <G( I )*G< I - 1 ) )*<DY/2. ) 

I F( IPRINT .NE. 0) 

SPRINT 600, I , Y ( I ) ,U( I ) , T ( I ) , F( I ) , SUMI , G( I ) , SUMG 
600 FORMAT(1X, I2.7F10.4) 

200 CONTINUE 

DELS = SUMI 
T H E T A = SUMG 

C CONTRIBUTION FROM WALL TO FIRST TUBE 

A = 0. 

B = Y(1)/DELTA 

CALL INCRS(MINF,REC, GAMMA, A, B, DELTA, L,EN,DELS1,THETA1, 

SI PR I NT , I ERR1 , I ERR2 ) 

I F ( I ERR 1 .LE. 2 .AND. IERR2 .LE. 2) GO TO 630 
PRINT 610, I ERR 1 , IERR2 

610 F0RMATI/10X, ‘ERROR IN SUBROUTINE CADRE WHILE CALCULATING*, 
SI OX , *C0NTR I BUT I ON FROM WALL TO FIRST TUBE:*, 

S10X, ‘ERROR NUMBERS :*,2I5) 

RETURN 

630 DELS = DELS + DELSl ‘(DELTA/L) 

THETA = THETA+ T H E T A 1 * ( D E L T A / L ) 

I F< IPRINT .NE. 0) 

SPRINT 660, DELSl, DELS , T H E T A 1 , T H E T A 
660 FORMAT(/10X,*DELS1 DELS THETA1 T H E T A* , 4 F 1 0 . 6 ) 

I FCDELTA .LE. Y(NBL)) RETURN 


C CONTRIBUTION FROM LAST TUBE TO DELTA 

A = Y(N8L)/L 
B - 1 .0 

CALL INCRS(MINF,REC, GAMMA, A, B, DELTA, L,EN,0ELS2, THETA 2, 
SIPRINT, IERR1 , IERR2) 

I F( IERR1 -LE. 2 .AND. IERR2 .LE. 2) GO TO 650 
PRINT 640, IERR1, IERR2 
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640 FORMAT(/10X, ‘ERROR IN SUBROUTINE CAORE WHILE CALCULATING*, 
S10X, ‘CONTRIBUTION FROM LAST TUBE TO DELTA:*, 

$1 OX , ‘ERROR NUMBERS :*,2I5) 

RETURN 

650 DELS * DELS ♦ DELS2 *(DELTA/L ) 

THETA « THETA+ T HE T A2* ( D E LTA/ L ) 

IF(I PRINT .NE. 0) 

SPRINT 670, DELS2 , THETA2 

670 FORMAT(/10X,*DELS2 T NET A2 : * , 2F 1 0 . 6 ) 

I F ( 0 E L T A .LE. Y < NBL ) ) RETURN 

C I OPT * 2: CALCULATE DELTASTAR AND THETA FROM Y«0 TO DELTA 

680 CONTINUE 
A * 0.0 
B * 1 .0 

CALL INCRS(MINF,REC, GAMMA, A, B, DELTA, L,EN, DELS , THETA , 

$ I PR I NT , I ERR 1 , IERR2) 

I F ( I E RR 1 .LE. 2 .AND. IERR2 .LE. 2) GO TO 700 
PRINT 690, IERR1,IERR2 

690 FORMATC/IOX, ‘ERROR IN SUBROUTINE CADRE WHILE CALCULATING*, 
S/10X,* FOR I OPT = 1 IN THE NAMELIST*, 

$/ 1 OX , ‘ERROR NUMBER: *.215) 

700 DELS = DELS • (DELTA/L) 

T H E T A = THETA* (DELTA/L) 

IF(! PRINT .NE. 0) PRINT 710, DELS, THETA 

710 F0RMAT(/10X,*DELS3 THETA3: *,2F10.6) 


RETURN 

END 
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SUBROUTINE I N CR S ( M I N F , R E C , G AHHA , A , B , D E L T A , l . E N , DE L S 1 , THETA 1 , 
S1PRI NT, IERR1 , IERR2) 

REAL MINF.L 
DIMENSION EPS ( 3 ) 

COMMON / AA/ C,ENP 
EXTERNAL FDEL, FTHETA 


GAMMA 


DELTA 


FREE-STREAM MACH NUMBER 

RECOVERY FACTOR 

RATIO OF SPECIFIC HEATS 

LOWER LIMIT OF INTEGRATION 

UPPER LIMIT OF INTEGRATION 

INDEX IN POWER LAW ( U/UE ) = ( Y /D ) ** ( 1 / N ) 

LENGTH SCALE 

BOUNDARY- LAYER THICKNESS (BY POWER-LAW FIT) 


DELS1 

THETA1 

IERR1 

IERR2 


VALUE OF DELSTAR/L FROM WALL TO NEAREST TUBE 
VALUE OF MOMENTM/L FROM WALL TO NEAREST TUBE 
NOT - 0, ERROR INTEGRATING DISPLACEMENT THICHNESS 
NOT = 0, ERROR INTEGRATING MOMENTUM THICKNESS 


C = (REC * (GAMMA - 1 . ) / 2 . ) * (MINF*MINF) 

ENP = EN 

IF(I PRINT .NE. 0) PRINT 100, M I N F , R E C , G AMMA , A , B , E N P , L , D E L T A 
100 FORMAT(8F1 0.6) 

E P S ( 1 ) = 1 . OE - 10 

E P S ( 2 ) = 1 . OE - 08 

ITEXT = 0 

CALL CADRE(A, B, FDEL, EPS, ITEXT, SUM, IERR) 

DELS1 = SUM 

IERR1 = IERR 

IFOPRINT .NE. 0) PRINT 200, SUM, IERR 
200 FORMATC10X.F10.6, 15) 


CALL CADRE(A,B, FTHETA, EPS, ITEXT, SUM, IERR) 
THETA1 = SUM 
IERR2 = IERR 

IF(IPRINT .NE. 0) PRINT 200, SUM, IERR 


RETURN 


s 


FUNCTION FDEL(YO) 
COMMON /AA/ C,ENP 


I F ( YD .NE. 

0 . ) U 

= YD**<1./ENP) 

IF(YD .eq. 

0. ) U 

= 0. 

T = 1 . 

+ C* ( 1 . 

- U*U) 

FDEL * ( 1 . 

- U/T) 


PRINT 100, 

YD , U , T , 

FDEL.C.ENP 

FORHATt 10X 

, 6 F 1 0 . 6 ) 



RETURN 

END 

FUNCT ION FTHETA(YD ) 

COMMON /AA/ C , E N P 

I F(YD .ME . 0. ) U * Y D * * ( 1 ,/ENP> 

I F ( Y D .EO. 0. ) U =0. 

T ■ 1 . ♦ C* < 1 . • U*U) 

F T H E T A = <U/T)*(1.-U> 

RETURN 

END 

FUNCTION COMP(P,T) 

C CALCULATES COMPRESSIBILITY FACTOR Z FOR 

C PRESSURE P(PSIA) AND TEMPERATURE T(DEG K ) 

P * P / 1 4 . 7 
BO* 1.37 
B 1 = - 8 . 773 E • 0 2 
B2 = 4 . 70 3 E - 04 
B 3 = -1.386E-06 
B 4 = 1 .462 E - 09 
C 0 = 5.521 
Cl* - 1 . 986 E - 0 1 
C 2 = 7.817 E ■ 04 
C3= -1.258E-06 
C 4 = 5 . 333 E - 1 0 

D = B0 + B1*T + B 2 * T * T +B3*(T**3.) +B4*(T**4.) 
E=C0+C1*T +C2*T*T +C3*(T**3.) +C4*(T**4.) 
2=1. -P* EXP(D) -<P*P)* EXP(E) 

C OM P = Z 
P = P * 1 4 . 7 
RETURN 
END 


GIVEN 
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by 



log(y/L) 


Figure 2: Determination of edge of boundary-layer and constant 

•n' in the power-law velocity profile (Uary 1 / 0 ) . 



Y/DELTR 



Figure 3 : Typical plot of measured boundary-layer velocities 

and power-law fit. 
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